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Abstract 

We consider a diffuse interface model for incompressible isothermal mixtures of 
two immiscible fluids with matched constant densities. This model consists of the 
Navier-Stokes system coupled with a convective nonlocal Cahn-Hilliard equation 
with non-constant mobility. We first prove the existence of a global weak solution 
in the case of non-degenerate mobilities and regular potentials of polynomial growth. 
Then we extend the result to degenerate mobilities and singular (e.g. logarithmic) 
potentials. In the latter case we also establish the existence of the global attrac- 
tor for the corresponding semiflow in dimension two. Using a similar technique, 
we show that there is a global attractor for the convective nonlocal Cahn-Hilliard 
equation with degenerate mobility and singular potential in dimension three. 
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1 Introduction 



Model H is a diffuse interface model for incompressible isothermal two-phase flows which 
consists of the Navier-Stokes equations for the (averaged) velocity u nonlinearly coupled 
with a convective Cahn-Hilliard equation for the (relative) concentration difference <p (cf., 
for instance, [6], E2, EH EDI EDI El])- The resulting evolution system has been studied by 
several authors (see, e.g., [U El QUI ESI ESI EDI ESI 1121 SI]) and references therein, cf. also 
0, ED ED] for models with shear dependent viscosity). In the case of matched densities, 
setting the constant density equal to one, the system can be written as follows 

u t - isAu + (u ■ V)w + Vtt = /iVy? + h infix(0,T) (1.1) 
div(w)=0 infix(0,T) (1.2) 
(ft + u ■ Vy? = div(m(<p)V//) infix(0,T) (1.3) 

H = -aAip + -F'(<p) in Q x (0, T), (1.4) 
a 

where Q C M d , d = 2, 3, is a bounded domain. Here z/ > is the viscosity (supposed to 
be constant for simplicity), ir is the pressure, h is a given (non-gradient) external force, 
m is the so-called mobility and a > is related to the (diffuse) interface thickness. 

A more realistic version of the Cahn-Hilliard equation is characterized by a (spa- 
tially) nonlocal free energy. The physical relevance of nonlocal interactions was already 
pointed out in the pioneering paper jlD] (see also [HI 4.2] and references therein). Though 
isothermal and nonisothermal models containing nonlocal terms have only recently been 
studied from the analytical viewpoint (cf., e.g., E21 [22J [231 EH ESI [33] and their ref- 
erences). The difference between local and nonlocal models consists in the choice of the 
interaction potential. The nonlocal contribution to the free energy has typically the form 
J n J(x,y) \<fi(x) — (f(y)\ 2 dy with a given symmetric kernel J defined on Q x Q; its local 
Ginzburg-Landau counterpart has the form (a/2)\'Vip(x)\ 2 with a positive parameter a. 
The latter can be obtained as a formal limit as m — > oo from the nonlocal one with the 
choice J(x,y) = m d+2 J(\m(x — y)\ 2 ), where J is a nonnegative function with support in 
[0, 1]. This follows from the formula formally deduced in [35] 

2 

dz 



m d+2 J(\m(x-y)\ 2 )\ V (x)- V (y)\ 2 dy= / J(|z| 2 ) 



j_ 

m 



J{\z\ 2 ) (V^(x),^) 2 dz = -\V<p{x) 

for a sufficiently regular <p, where a = 2/d J Rd J(\z\ 2 )\z\ 2 dz and Q m (x) = m(Q — x). 
Here we have used the identity J Rd J(\z\ 2 ) (e, z) 2 dz = 1/d L d J(\z\ 2 )\z\ 2 dz for every unit 
vector e G M. d . As a consequence, the Cahn-Hilliard equation fll.3p - fll.4l) can be viewed as 
an approximation of the nonlocal one. 
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Nonlocal interactions have been taken into account in a series of recent papers (see 
[TTj IT5| HH [T7]) where a modification of the model H with matched densities has been con- 
sidered and analyzed . More precisely, a system of the following form has been considered 

(cf. nu) 

<p t + u ■ Vp = div(m(<p)V//) infix(0,T) (1.5) 

fi= atp- J*ip + F'(ip) infix(0,T) (1.6) 

u t - uAu + (u ■ V)u + Vtt = fiVf + h infix(0,T) (1.7) 

div(u)=0 infix(0,T) (1.8) 

?E = o, u = 0, on90x (0, T) (1.9) 
an 

ii(0) = uo, </?(0) = <po, in f2, (1-10) 

where n stands for the outward normal to dQ, while Uq and po are given initial conditions. 
Here the interaction kernel J : lR d — > M is an even function and a(x) = / J(x — y)dy. 

Nonlocal system (jl.5B - fll.8p is more challenging with respect to (ll.ll) - (ll.4p . even in 
dimension two. One of the reasons is that ip has a poorer regularity and this influences 
the treatment of the Navier-Stokes system through the so-called Korteweg force uV</? (see, 
for instance, [TTI Remark 8]). Due to this difficulty, only the constant mobility case has 
been considered so far (though viscosity depending on tp has been handled). On the other 
hand, in the rigorous derivation of the nonlocal Cahn-Hilliard equation done in [23] the 
mobility depends on tp and degenerates at the pure phases (for the local Cahn-Hilliard 
equation see [131 [39] an d references therein). 

In this contribution we want to generalize some of the existing results on the so-called 
nonlocal Cahn-Hilliard-Navier-Stokes system to the case of nonconstant mobility. As we 
shall see, such an extension requires extra efforts which are not merely technical. 

In the seminal paper [13] the authors established the existence of a weak solution 
to the local Cahn-Hilliard equation with degenerate mobility and singular potentials en- 
dowed with no-flux boundary conditions. This result was then extended to the standard 
Cahn-Hilliard-Navier-Stokes system in [10]. The nonlocal Cahn-Hilliard equation with 
degenerate mobility and logarithmic potential was rigorously justified and analyzed in 
[23] (see also [25] and references therein). In particular, in the case of periodic boundary 
conditions, an existence and uniqueness result was proven in [23]. Then a more gen- 
eral case was considered in [22J . More recently, the convergence to single equilibria was 
studied in [36l EZ] (cf. also [21] for further results). 

Here we generalize the strategy devised in [13] to the nonlocal case by first taking a 
non-degenerate mobility m and a regular potential F with polynomial growth. We prove 
the existence of a global weak solution which satisfies an energy inequality (equality if 
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d — 2). This result extends [15] and allows us to construct a rigorous approximation of 
the case where m is degenerate and F is singular (e.g. logarithmic). Therefore we can 
pass to the limit and obtain a similar result for the latter case. In addition, since the 
energy identity holds in two dimensions, we can construct a semiflow which possesses a 
global attractor by using Ball's method (see [7]). By means of the same approach we also 
show that the convective nonlocal Cahn-Hilliard equation with degenerate mobility and 
singular potential has a unique solution and it posseses a global attractor (even if d = 3). 
Note that this result entails, in particular, that the nonlocal Cahn-Hilliard equation which 
has been obtained as hydrodynamic limit in [23] possesses a global attractor. We point 
out that uniqueness of solutions is still an open issue in the local case. 

Let us notice here that the main difficulty encountered while dealing with the degener- 
ate mobility case is that the gradient of the chemical potential /i in (11.61) can no longer be 
controlled in any L p space. Hence, in order to get an existence result a suitable notion of 
weak solution need to be introduced. More precisely, in this new formulation the gradient 
of \x does not appear anymore (cf. Definition [2] in Section m see also [T3j). It worth 
observing that, in the present case, our main Thm. [2] does not require the (conserved) 
mean value of the order parameter ip to be strictly in between —1 and 1, but | L ip \ < \Q\ 
suffices. Thus the model allows pure phase solutions for all t > 0. This was not possible 
in the case of constant or strongly degenerate mobility (cf. Remarks [6] and [7J for further 
comments on this topic). 

We conclude by observing that it would be particularly interesting (and challenging) to 
extend the present and previous results on nonlocal Cahn-Hilliard-Navier-Stokes systems 
to the case of unmatched densities (for the local model see [3j H] and references therein) 
or to the compressible case (cf. [5]). 

The plan of the paper goes as follows. In Section [2] we introduce the notation and the 
functional setting. Section [3] is devoted to prove the existence of a global weak solution 
satisfying a suitable energy inequality (equality in the 2D case) when the mobility is 
non-degenerate and F is regular. In Section HI using a convenient approximation scheme, 
we extend the proven result to the case of degenerate mobility and a singular F (e.g. 
of logarithmic type). Some regularity issues for ip and fi are also discussed. Section [5] 
is devoted to the existence of global attractor in the two dimensional case. Finally, in 
Section [61 we consider the convective nonlocal Cahn-Hilliard equation with degenerate 
mobility. We deduce the existence of a global weak solution from the previous result for 
the coupled system. Then, even in dimension three, we establish the uniqueness of weak 
solutions as well as the existence of a (connected) global attractor. 
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2 Notation and functional setting 



We set H := L 2 (Q) and V := if 1 (£7), where fl is supposed to have a sufficiently smooth 
boundary (say, e.g., of class C 1,1 ). If X is a (real) Banach space, X' will denote its dual. 
For every / G V we denote by / the average of / over Q, i.e., / := 1). Here \Q\ 

is the Lebesgue measure of Q. Let us introduce also the spaces Vq := {v G V : v = 0}, 
V ' := {/ e V : / = 0} and the operator A : V V , A e C(V, V) defined by 



(Au,v):= / Vm-Vw Vu,?;gV. 
in 



We recall that A maps V onto V ' and the restriction of A to Vo maps Vq onto V^' isomor- 
phically. Let us denote by Af : Vq — > V the inverse map defined by 

ANf = f, V/ G V ' and MAu = u, Vm G Vq. 

As is well known, for every / G Vq', M f is the unique solution with zero mean value of 
the Neumann problem 

—Am = /, in Q 
g = 0, on dQ. 

Furthermore, the following relations hold 

(Au,Aff) = (f,u), WueV, v/ey; 

(f,Ng) = (gMf) = [ v(Aff) ■ v(M g ), Vf,g g v '. 

Ju 

Recall that A can be also viewed as an unbounded operator A : D(A) C H — )■ if from 
the domain -D(y4) = {0 G H 2 (Vt) : dcfi/dn = on <9f2} onto if and that A/" can also be 
viewed as a self-adjoint compact operator M = A" 1 : H H in H. Hence, the fractional 
powers A r and J\f s , for r, s > 0, can be defined through classical spectral theory. 

We shall repeatedly need the standard Hilbert spaces for the Navier-Stokes equations 
with no-slip boundary condition (see, e.g., [13] ) 

-L 2 (n) d 



G div := {u G C^{Sl) d : div(u) = 0} \ V div := {u G H^O)* : div(ti) = 0}. 

We denote by || • || and (•, •) the norm and the scalar product on both H and Gdiv, 
respectively. Instead, Vd% v is endowed with the scalar product 

(u, v) Vdiv = (Vm, Vv), Wu, v G V div . 

In the proof of Theorem [1] we shall introduce the family of the eigenf unctions of the 
Stokes operator S with no-slip boundary condition. We recall that S : D(S) C G d i V ~^ 
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G div is defined as S := -PA with domain D{S) = H 2 {Q) d nV div , where P : L 2 {Q) d -»■ 
is the Leray projector. Notice that we have 

(Su,v) = (u,v) Vdiv = (Vu, Vi>), WueD(S), WveV div , 

and : Cr^u — >■ G^iu is a self-adjoint compact operator in G^. Thus, according with 
classical results, S possesses a sequence of eigenvalues {Xj} with < Ai < A 2 < • ■ ■ and 
Xj — > oo, and a family {wj} C D(S) of eigenfunctions which is orthonormal in Gdiv Let 
us also recall the Poincare inequality 

Ai||w|| 2 < ||Vm|| 2 , Vn G V div . 

The trilinear form 6 which appears in the weak formulation of the Navier-Stokes equa- 
tions is defined as follows 

b(u,v,w) = / (u • V)u • w, \fu,v,wEVdiv 
Jq 

We recall that we have 

b(u,w,v) = —b(u,v,w), \/u,v,w EVdiv 

For the basic estimates satisfied by the trilinear form b the reader is referred to, e.g., |13] . 

Finally, if X is a (real) Banach space and r 6 R, we shall denote by L^ b (r, oo; X), 
1 < p < oo, the space of functions / G Lf oc ([r, oo); X) that are translation bounded in 
Lf oc ([r, oo); X), i.e. such that 

rt+l 

ll/ll^(.oo.x) : = su P / \\f(s)\\ P x ds <oo. (2.1) 

3 Non-degenerate mobility 

Let us first consider the case where the mobility m does not degenerate, i.e. m satisfies 
the following assumption 

(HI) m G (7 J J C (R) and there exist mi,m 2 > such that 

mi < m(s) < m 2 , Vs G R. 

The other assumptions we need are the ones on the kernel J, on the potential F and the 
forcing term h which are the same as in [11] 



6 



(H2) J(- — x) G W 1,l (Q) for almost any ifll and satisfies 

J(x) = J(—x), a{x) := / J(x — y)dy > 0, a.e. iGO, 



a* := sup / | J(x — < oo, 6 := sup / | VJ(x — y)\dy < oo. 

xen Jn xen Jn 



(H3) F G Cf^(M) and there exists cq > such that 

F"0) + o(z) > c , Vs G R, a.e. xGil. 

(H4) There exist C\ > (a* — a*)/2, where 

a* := inf / J(x — y)dy, 
xenj n 

and C2 G M such that 

F(s) > ds 2 - c 2 , Vs G E. 
(H5) There exist C3 > 0, C4 > and r G (1, 2] such that 

\F'(s)\ r <c 3 \F(s)\+c 4 , VsGl. 

(H6) h G L 2 (0, T; V^J for all T > 0. 

Some further regularity properties of the weak solution can be established by using the 
following assumption 

(H7) F G C?£(R) and there exist c 5 > 0, c 6 > and p > 2 such that 

F"(s) + a(x) > c 5 |s| p " 2 - c 6 , VsGl, a.e. iGll. 

Remark 1. A well known example of potential F and kernel J which satisfy the above 
conditions is given by F(s) = (s 2 — l) 2 and J(x) = js\x\ 2 ~ d , if d = 3, and J(x) = 
—ji log |x|, if d = 2, where j'2 an d J3 are positive constants. 

Before stating the main result of this section, let us recall the definition of weak 
solution to system (jl.5j) — (jl.lOj) . 

Definition 1. Let uq G Gdiv, <Po G H such that F(ipo) G L 1 (Q), and < T < 00 be given. 
Then, a couple [u,<p] is a weak solution to fll.5p - fll.10l) on [0,T] if 

u G L°°(0, T; G dlv ) n L 2 (0, T; t^), <^ G L°°(0, T; iZ) n L 2 (0, T; V), (3.1) 
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u t eL^(0,T;V^ lv ), ^6L 4/3 (0,T;y'), ifd = 3, (3.2) 
Ui G^(0,T;4), p t eL 2 - s (0,T;V') 7,8 e (0,1), ifd=2, (3.3) 
/j:=o^-J^ + F^)Gl 2 (0,T;K), (3.4) 

and the following variational formulation is satisfied for almost any t G (0, T) 

(<p t , </>> + (m(<p)Vn, W>) = (ucp, VV), V0 G V, (3.5) 
(u t ,i;) + i^(V«, Vu) + b(u,u,v) = -(ipVfj,,v) + (h,v), VveV div , (3.6) 

together with the initial conditions fll.lOp . 

Remark 2. It is easy to see that u G C w ([0, T]; G^™) and G C w ([0, T]; if). Hence, the 
initial conditions (11.101) make sense. 

Theorem 1. Let u G Gdi V , y?o £ H such that F(ip ) G L 1 (f2) 7 and suppose that (Hl)- 
(H6) are satisfied. Then, for every given T > 0, there exists a weak solution [u,ip] to 
fll.Sp — fll.lOp in the sense of Definition^ satisfying the energy inequality 

S(u(t),<p(t)) + J (z/||Vu|| 2 + || v/m(^)V/i|| 2 )dr < £(u , y? ) + J (h(r),u)dr, (3.7) 

for every t > 0, where we have set 

£{u{t)^{t)) = \\\u{t)\\ 2 + l - [ [ J(x-y)( V (x,t)- V (y,t)) 2 dxdy + [ F(<p(t)). (3.8) 
A 4 Jn Jn Jn 

Furthermore, assume that assumption (H4) is replaced by (HI). Then, for every T > 
there exists a weak solution [u,(p] to fll.5p - fll.10l) on [0, T] corresponding to [uo,<po] (in 
the sense of Definition^) satisfying (13.11) . (13.21) and also 

ip G L°°(0, T; L p (f2)), (3.9) 
(fit G L 2 (0,T;V), if d = 2 or (d = 3andp>3), (3.10) 
u t EL 2 (0,T-X lv ), if d = 2. (3.11) 

Finally, assume that d = 2 and (H4) is replaced by (HI). Then, 
1. any weak solution satisfies the energy identity 

jfiu, ip) + v\\ Vwf + || ^m((p)Vn\\ 2 = (h(t),u), t > 0. (3.12) 
In particular we have u G C([0, oo); G^w), f G C([0, oo); H) and f n F(<p) G C([0, oo)). 
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2. If in addition h G -£4,(0, oo; V^ iv ), then any weak solution satisfies also the dissipative 
estimate 

£(u(t),<p(t)) < £(u ,(p )e- kt + F(m )\tt\ + K, Vt > 0, (3.13) 

where m = (ifio, 1) and k, K are two positive constants which are independent of 
the initial data, with K depending on Q, v, J, F and \\h\\ L 2 r 0oo -v' )■ 

tb^ ' ' div' 

Proof. The argument follows the lines of [HI Proof of Theorem 1] and is based on a 
Faedo-Galerkin approximation scheme. For the reader's convenience, we give a sketch of 
it. Let us first assume that <^ G D(B), where the operator B is given by B = — A + / 
with homogeneous Neumann boundary condition. We introduce the family {wj}j>i of the 
eigenf unctions of the Stokes operator S as a Galerkin base in V div and the family {ipj}j>i of 
the eigenf unctions of B as a Galerkin base in V. We define the n— dimensional subspaces 
W„ := (wi, • ■ ■ , w n ) and ^ n := (ipi, ■ ■ ■ , ip n ) and consider the orthogonal projectors on 
these subspaces in G div and H, respectively, i.e., P n := Py\> n and P n := P^, n . We then look 
for three functions of the form 

n n n 

«n(*)= E^H, ^(t)=^/t } (t)^, /in(t) = ^7i n) W^ 
k=l k=l k=l 

which solve the following approximating problem 

(<f' n , ip) + (m(v9„)V/i n , V-0) = (u n (p n , Vi/)) (3.14) 
(u' n , v) + v(Vu n , Vu) + b(Un, u n , v) = -(^V/i„, v) + (h n , v) (3.15) 
fj, n = P n (aip n - J*ip n + F'(ip n )) (3.16) 

^n(O) = <^0n, U n (0) = U 0n , (3.17) 

for every ip G \l/ n and every v G Wn, where ipo n = P n ^>o and wo™ = P n uo (primes denote 
derivatives with respect to time). In (13.15P {h n } is a sequence in C°([0, T}; Gdiv) such that 
h n — > h in L 2 (0,T; V^ itJ ). It is easy to see that this approximating problem is equivalent 
to solve a Cauchy problem for a system of ODEs in the 2n unknowns a 4 - n \ @!- (7}™^ 
can be deduced from (I3.16P ). Since F" and m are locally Lipschitz on R, the Cauchy- 
Lipschitz theorem ensures that there exists T* G (0, +00] such that system f l3.14p - fl3.17p 
has a unique maximal solution := (a[ n) , • ■ ■ ,ai n) ), b (n) := (/?J n) , • • ■ , 0) on [0,T*) 
with aW, b (n) G C^T^IT). 

By taking ip = fi n and t> = u n in f)3.14p and f l3.15p . respectively, and adding the 
resulting identities together, we get 

^-£(u n ,ip n ) + z/||Vw n || 2 + || v / "i(v5n)V/i n || 2 = (h n ,u n ), 
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where £ is defined as in (13. 8p . Integrating this identity between and t 

£(u n (t),(p n (t)) + J (v\\Vu n \\ 2 + ||A/m(v?n)VyU„|| 2 jdr = £{uon,(p 0n ) + J (h n (r),u n )dT. 

(3.18) 

Observe now that 

V 1 ■— "27 , 1 / NT 112 



(h n {j),u n )dT < - I \\Vu n \\ dr + 7^ / \\ h ™\\v^ v dr - 

J «/ 

On the other hand, taking (H2) and (H4) into account, we get 

111 f 

£(u n ,<p n ) = -\\u n \\ 2 + -\\Vaip n \\ 2 - -(ip n , J* if n ) + / F(ip n ) 

— — — J Q 

>^\Wn\\ 2 + ^ I aip 2 n - ];\\<Pn\\\\J * <Pn\\ + f (cifl - Pa) 

> ^ll^nll 2 + c'lll^nll 2 - c 2> 

where c' x = (a* — a*)/2 + c\ > and c' 2 = C2|f2|. Using also the convergence assumption 
for {u 0n }, {h n }, the fact that (po n —> ipo in H 2 (Q) (since ipo G D(B)) and the lower 
bound mi > for the mobility m (cf. (HI)), we first deduce that T* = +oo for every 
n > 1 (notice that |a^(t)| = ||w n (t)|| and |b^(t)| = H^nWH) and furthermore we get the 
following estimates which hold for any given < T < +oo 

\\u n \\L°°(0,T;G div )nL2(0,T;V div ) < C, (3.19) 
||<Ai||l/x>(0,T;ff) < C, (3.20) 
\\F(<p n )\\ L o° {0 ,T;Li m <C, (3.21) 

||V/i n ||L2( 0iT .#) < C, (3.22) 

with C independent of n. Henceforth we shall denote by C a positive constant which 
depends at most on ||wo||, \\fo\\, JqF^q), ||/i||i2(o,T;V'. ) an d on J, F, u, mi, Q, but they 
are independent of n. Instead, c will stand for a generic positive constant depending on 
the parameters of the problem only, i.e. on J, F, v, m\ and Q, but it is independent of 
n. The values of both C and c may possibly vary even within the same line. 
Let us now recall the estimate 

c 2 

||V/x n || 2 > ^||V^ n || 2 -c|K|| 2 , (3.23) 

which can be deduced as in [TTj Proof of Theorem 1] by multiplying fi n by — Aip n in L 2 
and integrating by parts. By means of (I3.23j) . from ( I3.20p and (I3.22p we get 

hn\\ L *(0,T;V) < C. (3.24) 
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Due to (H5), which in particular implies that < c|F(s)| + c, we have 

|/Z n |= [ F'((p n ) <c\\F(<p n )\\ LHn) + c<C, 
Jq 

and therefore, thanks to (I3.2ip . (I3.22p and to the Poincare-Wirtinger inequality, we deduce 

||/^n||L 2 (0,T;V) < C. (3.25) 

In addition, (H5) and (I3.2ip yield the following control 

\\F'(tp n )\\ L °°(0,T;Lr m <C. (3.26) 

Let us now derive the estimates for the sequences of time derivatives {u' n } and {<p' n }. 
By taking ipj := P n ip, foripeV arbitrary, as test function in ( I3.14p we obtain 

Assume first d = 3. We have 

Wn^)\ < "i2||V/i n ||||V^/|| + ||tt„||l,6 (n )3||<p n || L 3 ( n)||VV>j|| 

< c(||V/i n || + ||Vw„||||^ 

"II ^ 1 1 Ili6(r2) livvil 

< C(\\Vfin\\ + ||V« n ||||^ n ||v /2 )||VV||, W» e V (3.27) 

where in the last estimate we have used f)3.20p . Then, by f)3.19p . (I3.22p and (13.241) from 
(13T27j) we get 

Wn\\Li/'i(0,T;V) < C. (3.28) 

For d = 2, by arguing as above it is not difficult to infer the following bound 

||^nlk 2 - a (o,T;y) < C. (3.29) 

As far as the sequence {u' n } is concerned, by arguing exactly as in [TH Proof of 
Theorem 1] we can deduce the bounds 

\WJ\l*/3(o,T;V^ v ) < C> if d = 3, (3.30) 

\WJl^(o,T;V^) <C, if d = 2. (3.31) 

From fl3TT9l - fl3T20l . fl3T24D - fl3T26D and fl3T28l - fl3T3ll . using compactness results, we obtain 
for a not relabeled subsequence 

weakly* in L°°(0, T; G div ), weakly in L 2 (0, T; V div ), (3.32) 

u n — > u strongly in L 2 (0, T; Gdi V ), a.e. in Q T , (3.33) 
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u' n ^u t weakly inL 4 / 3 (0,T;\/j.J, d = 3, (3.34) 

u' n -± u t weakly in L 2 " 7 (0, T; V^J, d = 2, (3.35) 

y> n ->> 9? weakly* in L°°(0, T; if), weakly in L 2 (0, T; V), (3.36) 

ip n — >■ <y9 strongly in L 2 (0, T; if ), a.e. in Qr, (3.37) 

^ -> ^ weakly in L 2 ~ 5 (0, T; V"), (3.38) 

F'(yy) ->> F* weakly* in L°°(0, T; L p (fi)), (3.39) 

fin^fi weakly in L 2 (0, T; V). (3.40) 



where := fix (0, T). The pointwise convergence f)3.37p . the weak * convergence f )3.39p 
and the continuity of F' yield F* = F'(<p). By means of f l3.32p -f l3.40p we can now pass 
to the limit in the approximate problem f)3.14p - f)3.17p as done in [TU Proof of Theorem 
1] and obtain the variational formulation of system f ll.5p - fll.10p for functions u, <p and \x. 
In particular notice that, fixing n arbitrarily, for every \ £ C£°(0,T) and every ip 6 ty n 
we have 

/ (m(<f k )Vnk, Vi/j)xdr ->■ / (m(ip) V/i, S7ip)xdr, as fc ->■ oo, 
as a consequence of the weak convergence f)3.40p and of the convergence 

m(<£>ft) — > m(ip) strongly in L v (Qt), for all r£[l,oo), (3-41) 

ensured by (I3.37P and by Lebesgue's theorem. 

As far as the energy inequality is concerned, let us observe that (I3.40p and (I3.4ip imply 
that 

a/ 'm(ip n )V /i n ^m(if)Vfi weakly in L S (Q T ), for all s G [1,2). (3.42) 
But ( EHD yields 

||A/m(v9 n )V/i„|| L 2 ( Q T) < C, 
and hence f)3.42p holds also for s = 2. Therefore 

/ || ^/m(ip)'V fi\\ 2 dr < liminf / || a/ m(<y9 n )V/i„|| 2 dr. 
Jo n ^°° Jo 

This fact and the above convergences allow us to pass to the inferior limit in (I3.18P and 
deduce (13.71) . 

Let us now suppose that assumption (H4) is replaced by (H7). Then, the energy of 
the approximate solution can be controlled from below by 

£(u n ,(p n ) > c(||?I n || 2 + H^nH^n)) - c. 
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From (I3.18P we can then improve ( I3.20p . that is 



\\ i Pn\\L° o (0,T;LP(n)) < C, (3.43) 

which yields ( 13. 9p . The regularity properties (I3.10p and (13. lip follow by a comparison 
argument exactly as in [TTj Proof of Corollary 1] with reference to the weak formulation 
( 13. 5p . ( 13. 6 p and taking the improved regularity for yj and (HI) into account. 

Furthermore, if d — 2 and (H4) is replaced by (H7), we can take fi and u as test 
functions in the variational formulation (13. 5p . (13. 6p and argue as in [TTj Proof of Corollary 
2] in order to deduce the energy identity (I3.12p . As far as the the dissipative estimate 
( I3.13P is concerned, this can be deduced without difficulties by adapting the argument of 
[TTj Proof of Corollary 2] and taking (HI) into account. 

Finally, if (po £ H with -F(y?o) £ L 1 (f2), we can approximate (po with y?om G D(B) 
given by ip 0m := (J + B /m)^ 1 ^®. This sequence satisfies y?om — > <p in H and by exploiting 
assumption (H3) and the argument of [TT| Proof of Theorem 1] we can easily recover 
the existence of a weak solution, the energy inequality (13. 7p and, for d — 2, the energy 
identity ( I3.12p as well as the dissipative estimate ( I3.13p . □ 

It may be interesting to observe that another energy identity for d = 2 (inequality for 
d — 3) is satisfied by the weak solution of Theorem [TJ This energy identity will turn to 
be useful especially in the degenerate case in order to establish the existence of the global 
attractor in 2D. 

Corollary 1. Let the assumptions of Theorem^ be satisfied with (H4) replaced by (HI). 
Then, if d = 2 the weak solution z = [u, <p\ constructed in Theorem^ satisfies the following 
energy identity 

ljA\\ u \\ 2 + M 2 )+ I m( V )F'\ V )\W V \ 2 + f am(^)|Vyf + u\\Vu\\ 2 

= / m(<y?)(VJ* y? - y?Va) ■ Vy? + (cup - J * <p)u ■ Vy? + (h,u) , (3.44) 
Jn Jn 

for almost any t > 0. Furthermore, if d = 3 and if (H7) is satisfied with p > 3, the weak 
solution z satisfies the following energy inequality 



k\Ht)\\ 2 + h(t)\\ 2 ) + [ f m^)F"^)\W V \ 2 + [ fam(<p)\Vrf 
z Jo Jn Jo Jn 



+ v [ \\Vu\\ 2 < ^(|K|| 2 + ll^oll 2 ) + / [ m((p)(VJ *<p - <pVa) ■ Vy? 
Jo ^ Jo Jn 

+ (a<p- J*<p)u-V<p+ (h,u)dr, Wt > 0. (3.45) 

Jo Jn Jo 



13 



Proof. If d = 2, we can take ip = ip and v = u as test functions in the weak formulation 
(13. 5p and (13.61) . respectively. This is allowed due to the regularity properties (13.101) and 
(13. lip . Then, we add the resulting identities and we observe that, since by an Alikakos' 
iteration technique as in [HI Theorem 2.1] it can be shown that ip G L°°(Q) and since 
F G C 2 (M), we can write VF'(ip) = F"(ip)V<p. Then, on account of this identity and 
of (13.41) . the second term on the left hand side of (13.51) can be rewritten and (13.441) 
immediately follows. 

Let d = 3 and assume that (H7) holds with p > 3. Then, due to the second of (13.101) 
we can take ip = ip as test function in (13. 5p and get 

~j t y\\ 2 +(m{ip)Vv,Vip) =0. 

By integrating this identity between and t and by rewriting the term (m(ip) V/x, Vy?) as 
done above in the case d = 2 we obtain 

^)l| 2 + f [ m(ip)F'\ip)\Wip\ 2 + f [ am(<p)\V<p\ 2 
z Jo Jq. Jo Jn 

+ [ I m{ip)(ipVa- VJ*<p) ■ Vy = -||^o|| 2 . (3.46) 
Jo Jn 2 

On the other hand we take v = u n in the second equation (13.151) of the Faedo-Galerkin 
approximating problem and integrate the resulting identity to get 



^||u„(t)|| 2 + v j \\Vu n \\ 2 dr = l-\\u 0n \\ 2 - I (ip n Vfi n ,u n )dr + f (h,u n )dr 
Jo Jo Jo 

/ / [aip n - J * ip n )Vip n ■ u n + / (h,u n )dr. (3.47) 
Jo Jn Jo 
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2 II '"On || 2 



We can now pass to the limit in ( I3.47P and use the second weak convergence ( I3.32p and 
( I3.33p . In particular, observe that 

J J (aip n - J* if n )Vip n ■ U n = -J Q^P n Va - ifniVJ* ip n )^j ■ U n 

— > — j 0-ip 2 Va — ipiyj * ip) V u = j I (cap — J * ip)Vip ■ u, 
JQt ' Jo Jn 

where this last convergence is a consequence of the fact that, due to (13. 32 p . (13. 33 p . (I3.36j) . 
( I3.37P and interpolation, we have <p n — )■ ip strongly in L 3 (Q t ) and u n — >■ u strongly in 
L 3 {Qt) 3 (recall that Va G L°°). Hence we get 

\\\u(t)\\ 2 + u [ \\Vu\\ 2 < ^||m || 2 + I [ (cup- J*ip)u-Vip+ [ (h,u)dr. (3.48) 
2 io 2 Jq Jq J 

Summing (ERB]) and (E3HD we deduce $333]). □ 
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Remark 3. The theorems proven in this section still hold when the viscosity smoothly 
depends on (p (see [15J). Moreover, they allow us to generalize to the variable (non- 
degenerate) mobility the results obtained in [T5] on the longtime behavior. More precisely, 
the existence of the global attractor and the existence of a trajectory attractor in the cases 
d = 2 and d = 3, respectively. 

4 Degenerate mobility 

In this section we consider a mobility m which degenerates at ±1 and we assume that the 
double-well potential F is singular (e.g. logarithmic like) and defined in (—1,1). More 
precisely, we assume that m G C 1 ([— 1, 1]), m > and that m(s) = if and only if s = — 1 
or s — 1. Furthermore, we suppose that m and F fulfill the condition 

(Al) F G C 2 (-l,l) and 

mF" G C([-l,l]). 

We point out that (Al) is a typical condition which arises in the Cahn-Hilliard equation 
with degenerate mobility (see [131 1231 121122]). 

As far as F is concerned we assume that it can be written in the following form 

F = F l + F 2 , 

where the singular component F\ and the regular component F 2 G C 2 ([— 1, 1]) satisfy the 
following assumptions. 

(A2) There exist a 2 > 4(a* — a* — b 2 ), where b 2 := minj-^i] F 2 , and e > such that 

F?(s) > a 2 , Ws G (-1, -1 + e ] U [1 - e , 1). 

(A3) There exists eo > such that F^ is non- decreasing in [1 — eo, 1) and non-increasing 
in (-1,-1 + e ]. 

(A4) There exists Co > such that 

F"(s) + a(x) > c , VsG(-l,l), a.a. x G Q. (4.1) 

The constants a* and a* are given in (H2) and (H4), respectively, and the assumption 
on the external force h is still (H6). 
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Remark 4. It is easy to see that (Al)-(A4) are satisfied in the physically relevant case 
where the mobility and the double-well potential are given by 

m{s) = h(l - s 2 ) } F{s) = -|s 2 + °-({l + s) log(l + s) + {l-s) log(l - s)), (4.2) 

where < 9 < 9 C . Indeed, setting F x {s) := (6/2)((l + s) log(l + s) + (1 - s)log(l - s)) 
and F2(s) = — {9 c /2)s 2 , then we have mF" = k\9 > and so (Al) is fulfilled. Moreover 
Fi satisfies also (A2) and (A3), while (A4) holds if and only if inf^a > 9 C — 9. Another 
example is given by taking 

m{s) = k{s){l - s 2 ) m , F{s) = -k 2 s 2 + F^s) 

where k G C 1 ([— 1, 1]) such that < k 3 < k(s) < k± for all s G [—1,1], and F\ is a 
C 2 (— 1, 1) convex function such that 

F{'(s) = l(s)(l-s 2 r m , V S G(-1,1), 

where m > 1 and / G C 1 ([— 1, 1]). 

Remark 5. Note that (A4), which is equivalent to the condition inf^a > — inf F", 
is more general than [THl (A6)]. More precisely, [THl (A6)] implies (A4), while (A4) implies 
[TBI (A6)] provided we have inf (-1,1) (F" + F!f) = inf (-1,1) F" + minj-i^] F![. For example, 
consider the following double-well potential 

F(s) = -^s 2 -^+ 6 -{(l + s) log(l +s) + {l-s) log(l -a)), 

where < 9 < 9 C and < 9 < 9 2 . Then, it easy to see that [TBJ (A6)] is satisfied iff 
inf^ a > 9 C + 9 2 — 9, while (A4) requires the weaker condition inf^ a > 9 C + 9 2 — 2y/99 2 . 

If the mobility degenerates we are no longer able to control the gradient of the chemical 
potential fi in some L p space. For this reason, and also in order to pass to the limit in 
the approximate problem considered in the proof of the existence result, we shall have to 
suitably reformulate the definition of weak solution in such a way that \x does not appear 
anymore (cf. [T31). 

Definition 2. Let Uq G G^v, V^o £ H with F{ip ) G L x (fi) and < T < +00 be given. A 
couple [u,<f] is a weak solution to f ll.5p - fll.10p on [0,T] corresponding to [uq, (fo) if 

• u, ip satisfy 

u G L°°(0, T; Gdiv) n L 2 (0, T; V dm ), (4.3) 
u t EL^(0,T-X lv ), if d = 3, (4.4) 
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u t eL 2 (0,T;V^ v ), if d = 2, (4.5) 

<p G L°°(0, T; H) fl L 2 (0, T; V), (4.6) 

Vt eL 2 {S) ) T-V') 1 (4.7) 

and 

tpeL°°(Q T ), \ip(x,t)\<l a.e. (x,t) G Q T := x (0,T); (4.8) 

• /or every ^ 6 V, every t> G and /or almost any t G (0, T) we /iave 

(y>t,V) + / rn{<p)F"(<p)V(p- Vip+ / m(^)aV^ • 
</n </q 

+ / m{<p){<pVa - VJ* y>) ■ = (up, W), (4.9) 
in 

(«t, u) + r/(Vu, Vr) + b(u, u, v) = ((aip -J* ip)Vp, v) + (h, v); (4.10) 

• the initial conditions u(0) = u$, <f(0) = fo hold (cf. Remark^). 

We now state the existence result for the degenerate mobility case. To this aim we 
need to introduce the function M G C 2 (— 1, 1) defined by 

m{s)M"{s) = 1, M(0) = M'(0) = 0. 

Theorem 2. Assume that (Al)-(A4) and (H2), (H6) are satisfied. Let uq G fo G 

L°°(Q) such that F(ipo) G L x (f2) and M((po) G L 1 (Q). Then, for every T > t/iere exists 
a weak solution z := to (1 1.5ft - (1 1.1 Oft on [0,T] corresponding to [u , (fo) in the sense 

of Definition^ such that Tp(t) = Tp^ for all t G [0, T] and 

peL oo (0,T;ZP(n)) l (4.11) 

where p < 6 /or d = 3 and 2 < p < oo /or d = 2. In addition, if d = 2, the weak 
solution z := [it, p] satisfies the energy equation (13.44j) . and, if d = 3, z satisfies the 
energy inequality (13. 45 p . 

Proof. Let us consider the approximate problem P e : find a weak solution [u e , (p e ] to 

+ u e ■ Vip e = div(m e (p e )V/x e ), (4.12) 
u' e - uAu e + (u t ■ V)u t + Vvr e = u e Vy? e + h, (4.13) 

He = We - J * fe + F' e (<f e ), (4.14) 

div(u e ) = 0, (4.15) 
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?tl = o, w e = 0, on dn, (4.16) 
an 

w e (0) = u , (fe(0) = <Po, in fi. (4.17) 

Problem P e is obtained from f 1 1 . 5 p - f TTTTD]) by replacing the singular potential F with a 
smooth potential F t and the degenerate mobility m with a non-degenerate one m e . In 
particular, F t is represented by In particular, F e is represented by 



F e = F u + F 



2c- 



where F u and F 2e are defined by 

Fl{s) = \ ./•;'(.). H<i-r (i.i.s) 




^OOH M<l-c ( 4 - 19 ) 



with Fie(O) = Fi(0), ^(0) = F x '(0), F 2e (0) = F 2 (0), i^(0) = F&0). 
The approximate non-degenerate mobility is given by 

{m(l — e), s > 1 — e 

|s|<l-e (4.20) 
m(-l + e), s < -1 + e. 

Assumption (A3) implies that 

F le (s) <F 1 (s), Vsg(-1,1), Vee(0,e o ]. (4.21) 
On the other hand, from the definition of F 2e we have, for s > 1 — e, 

F 2e (s) = F 2 (l - 6) + F 2 '(l - e) (5 - (1 - e)) + - 6) (a - (1 - e)) 2 , (4.22) 

(a similar expression holds for s < — 1 + e) and F 2e (s) = F 2 (s) for |s| < 1 — e. Since 
F 2 G C 2 ([— 1, 1]), then we deduce that there exist two positive constants L±, L 2 , which are 
independent of e, such that 

\F 2e (s)\ < L lS 2 + L 2 , Vs G E. (4.23) 
Therefore, by using the assumption on the initial datum <p and f)4.2ip . (I4.23[) we have 

[ F £ (<p ) < f F l { VG ) + L i y Q \\ 2 + L 2 < oo, VeG(0,e ]. (4.24) 
Jn Jn 

18 



Furthermore, we can see that by assumption (A2) there is a So G R such that 

F e (s) > Ss 2 - S , Vs Gf, Ve G (0, e ], (4.25) 

where < 5 < a 2 /8 + b 2 /2. Indeed, (A2) implies that there exist G R such that 

F{(s) > Ox for all s G [1 - e , 1), F{(s) < aj~ for all s G (-1, -1 + e ], and Fi(s) > a for 
all s G (—1,1). Hence, by using definition (I4.18p . for s > 1 — e we get 

F„( S ) = F,(l - e) + Fl(l - e){s - (1 - e)) + ifftl - e)( S - (1 - e)) 2 
>a i + |( S -(l- £ )) 2 >|^ + a ;-|, (4.26) 

provided that e G (0, e ]. For \s\ < 1— e we have F le (s) = Fi(s) > a > (a 2 /8)s 2 — a 2 /8+a - 
For s < — 1 + e we can argue as in ()4.26|) . On the other hand, from f)4.22p we have 

F 2t (s)>(S-^s 2 -S 1 , VsGR, VeG(0,l), (4.27) 

where S is taken as above and 8\ is a nonnegative constant (depending on b := minr_n] F 2 , 
6i := min [ _i i i ] and b 2 ). Combining (S2SD with (jOZD we deduce f!4~25]) . 
Moreover, (A4) immediately implies that, for e G (0, eo], there holds 

F e "(s) + a{x) > c , Vs G R, a.e. iGfl. (4.28) 

Indeed, recall that for s ^ ±1 =f e we have f£(s) = F"(±l =F e), F£(s) = i^'(±l =p e), and 
for |s| < 1 - e we have F{' € (s) = F?(s), F£(s) = Fjfts). 

We can now check that assumptions (H1)-(H6) of Theorem [1] for Problem P e (with 
a fixed e G (0, eo]) are satisfied. In particular, due to (A2), we can choose S such that 
(a* — a*)/2 < 5 < a 2 /8 + b 2 /2 and so (H4) is ensured by (14.251) . while (H5) is satisfied with 
r = 2. Theorem [1] therefore entails that Problem P e admits a weak solution z e = [w<:,<^<:] 
satisfying (I3.ip -( 13~U]) and the energy inequality 

S e (z e (t))+ J (v\\Vu e \\ 2 + || v/m e (^)V/i £ || 2 )dr < £ e (z ) + J (h(r),u e (r))dT, (4.29) 
for every t > 0, where 

111 /" 

^e(z e ) := -||mJ 2 + -Hva^f - -(^ 6 , J * Ve) + J F e (<p e ). 

Now, using ( I4.24p and ( I4.25p . from (I4.29P we immediately obtain the following uniform 
(with respect to e G (0,e ]) estimates 

\\Ue\\L°°{0,T;G div )nL2{0,T;V div ) < C, (4.30) 
H^elU^CO.T^) < C, (4.31) 
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\\\/m e (ip e )W ^ e \\ L 2 {0)T . H) < C, (4.32) 

where the positive constant C can possibly change from line to line but is always inde- 
pendent of e. Now, test (14.121) by ip = M' e ((p e ), where M t is a C 2 function such that 
m e (s)M'J(s) = 1 and M e (0) = M' e (0) = 0. We get 



d_ 

It 



[ M e {tp e ) + [ m e {<p e )Vfi e ■ M'JicpjVcpe = [ u e cp e -M' e '{cp e )Vcp e = 0, (4.33) 
Jn Jn Jn 



where the last identity is due to (14.151) . Therefore (I4.33P yields 

^ f M e (if e ) + f V/4 • Vcp e = 0. 



dt 

On account of (14.141) . we obtain 



n Jn 



j t ^ M e (<p e ) + J ((a + F e 'V)) | V^ e | 2 + ^ e Va • V<^ e - V J * y? e • Vy? e ) = 0. (4.34) 



By using (14. 28p and (I4.3ip . from (I4.34p we hence infer 

/ M e (^) + ^||V^|| 2 < CJ ||^|| 2 <C, 
Jn 1 



d_ 

dt 



where cj is a positive constant depending on J only. Therefore, on account of the fact 
that for e small enough we have M e (s) < M(s) for all s G (—1, 1), and recalling that 
J n M({p ) < oo, we deduce the bounds (see also (14.311) ) 

lke|U 2 (0,T;V0 < C, (4.35) 

||M e (^ e )|| ioo( o ir;X1( n)) < C. (4.36) 

Let us now establish the estimates for <p' e and u' e . Let us first consider the case d = 3 
and start from (p' e . From the variational formulation of (I4.12p we get 

1(^,-0)1 < ll^e(^e)V// e ||||V0|| + |MU 8 (0)|IV 7 e||ii 3 (n)||VV'|| 

< c||v/m e (^jV^J||V^|| + c||VM e ||||^ir / l^|li / 6 2 (n) ||V^|| 

< C(||Vm e (^)V/i e || + \\Vu e \\ye\\\l 2 W> G V, (4.37) 

where we have used (14.311) . Then, on account of (I4.30p . (I4.32p and (I4.35p . from (14.371) we 
obtain 

lkellL4/3( ,T;V') < (4.38) 

In order to deduce a bound for u' t , observe first that we have 

(// e V</?e, V) = ((aipe - J* V?e) Vy? e , v) , Wv G V div . 
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Thus, by f)4.3ip . we have that 



\(/j, e V(p e ,v)\ < \\aip e -J* ¥> e ||i3(n)||Vy> e ||||v||z,6( n ) 

< cj\\<p e \\ LH n)ye\\v\\v \\ Vdiv < cj||<^ e || 1/2 ||^||2 6 2 {n) ||v? e || y ||'y \\ Vdiv 

< C\\<p e \$ 2 \\v\\ Vdiv , (4.39) 

Hence, (I4.35P entails the estimate for the Korteweg force in (I4.13P 

ll/ i eVV2 e || L 4/3( 0T y'. )<C. 

V 1 d%v ' 

The estimates for the other terms in the variational formulation of (14. 13j) can be deduced 
by standard arguments for the 3D Navier-Stokes and finally we obtain 

KIIl4/3(o,t ; v^) < C. (4.40) 

In the case d = 2, by arguing as in (I4.37P and ( I4.39P and again using f l4.30p - fl4.32p and 
f!4.35p . it not difficult to see that the following estimates hold 

<C, V 7 ,5g(0,1). (4.41) 

From the estimates above and using standard compactness results we deduce that there 
exist u G L°°(0, T; G div ) n L 2 {0, T; V div ) and ip G L°°(0, T; H) n L 2 (0, T; V) such that, up 
to a (not relabeled) subsequence we have 

u e ^u weakly* inL°°(0,T; G div ), weakly in L 2 (0, T; V div ), (4.42) 

u e — » u strongly in L 2 (0, T; G div ), a.e. in Q T , (4.43) 

u e -± u t weakly in L 4/3 (0, T; V^ v ), if d = 3, (4.44) 

< - « t weakly in L 2 ^(0,T; V^), 7 € (0, 1), if d= 2, (4.45) 

<p e ->> up weakly* in L°°(0, T; if), weakly in L 2 (0, T; V), (4.46) 

ip e ^)-ip strongly in L 2 (0, T; H), a.e. in Qt, (4.47) 

<p' t -> ^ weakly in L 4/3 (0, T; 7'), if d = 3, (4.48) 

<p' e p t weakly in L^O, T; V'), 5 G (0, 1), if d = 2. (4.49) 

Let us now show that \ip\ < 1 almost everywhere in Qt- In order to do that we can argue 
as in [131 Proof of Theorem 1] (see also [TC)| Proof of Theorem 2.3]). More precisely, in 
Proof of Theorem 1] the following estimates are established 

"■<*> s ^hr^- 1)2 - Vs>1 < 



21 



Therefore 

<2m(l-e) / M e (<fi e ) + 2m(-l + e) / M e (^ e ) 

<2max(m(l-e),m(-l + e)) / M e (v? e ). 

Now, by using (14.361) . the fact that m(±l =)= e) — > as e — > and the generalized Lebesgue 
theorem, we obtain 

/ (|^| - i) 2 = for a.a. t G (0,T), 

so that 

\(p(x,t)\<l, a.e. in fi, a.a. te(0,T). (4.50) 

We now have to pass to the limit in the variational formulation of the approximate 
problem (14. 12j) (14. 17j) in order to show that [u, ip] is a weak solution to (jl.5j) — (ll.lOj) ac- 
cording to Definition [2J It is easy that the variational formulation of (I4.12p - (l4.17p gives 

&' e ,ip) + / m e {(p e )F?(<p e )V<p e -Vil>+ I m e {if e )aVif e -V^ 
Jn Jn 

+ / m e (<p e )(<p e Va - VJ * <fi e ) ■ Vlf) = (U e <fi e , Vip), (4.51) 
Jn 

(u' e , v) + u(Vu e , Vv) + b(u e , u e , v) = ({aifi e - J * ifi e )V<fi e , v) + (h, v), (4.52) 

for every if) G V, every v G Vdiv and for almost any t G (0, T). 

On account of (Al) we have \m(s)F"(s)\ < a for every s G [—1, 1] and for some positive 
constant a. This immediately implies that 

K(s)F e "(s)| < a, Vs G R, Ve G (0, 1). (4.53) 

Furthermore, by using the continuity of mF" on [—1, 1] (cf. (Al)) and (14.471) . it is not 
difficult to see that 

m e (ip t )F"(ip e ) -)> m(<f)F"((fi), a.e. in Q T . 

Therefore, we obtain 

m e {(p e )F^(p e )^m{(p)F"((p), strongly in L r (Q T ), VrG[2,oo). (4.54) 
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Notice also that (Ell, the continuous embedding L°°(0, T; L 2 (tt)) ni 2 (0,T;L 6 (O)) <^ 
L 10/3 (Q T ) and (Q7|) imply that, for d = 3, we have 



^ e ->• strongly in L"(Q T ), for all 2 < s < 10/3. (4.55) 

For d = 2, by the Gagliardo-Nirenberg inequality the continuous embedding L°°(0, T; L 2 (f2))fl 
L 2 (0,T; V) L 4 (Q T ) holds, thus we get 

V? e -)> strongly in L S (Q T ), for all 2 < s < 4. (4.56) 

Furthermore, by Lebesgue's theorem we also deduce 

m e {ip e ) — > m(<f) strongly in L s (Qt), for all 2 < s < oo. (4-57) 

We can now multipy (14.511) . (14.521) by x,cj G C£°(0, T), respectively, and integrate the 
resulting identities with respect to time from and T. The convergences (I4.42I) - (I4.49I) 
and (I4.54I) - (I4.57I) are enough to pass to the limit and to deduce that u and ip satisfy 
the variational formulation (14.91) and (14. 10p . Observe that in passing to the limit we can 
initially take ip in C 1 (Q) and then prove by density that (14. 9 p holds also for all ip in V. 

In order to prove that z := [u, ip] is a weak solution according to Definition [21 we are 
only left to show that (14. 5p and (14.71) hold. To this aim let us first see that ip satisfies 
(14. lip . Indeed, this is a consequence of f)4.5Up and of the fact that, since ip G L 2 (0,T; V), 
ip is measurable with values in L 6 (Q) when d = 3, and measurable with values in L P (Q), 
for all p G [2, oo), when d = 2. Recalling the improved regularity (14. lip , and also on 
account of (I4.50p . it is now immediate to get (14. 5 p and (14. 7p by a comparison argument 
in the weak formulation ( I4.9p . ( I4.10p . In particular, notice that the contribution from the 
convective term in ( 14 .9p can now be estimated as follows 

\{ucp, w)l < IMIz^IMIlWIWH < C||V«||||V^||, vv> g V. 

Finally, if d — 2, the energy identity (13.44j) can be obtained by taking i\) = ip and 
v = u (see (14. 7p and ( I4.5P ) as test functions in the weak formulation ( 14. 9 p and (I4.10p . 
respectively, and adding the resulting identities. In the case d — 3, we choose first ip = ip 
as test function in (14. 9 p (cf. (14. 7p ) and integrate the resulting identity to get (13. 46 p . Then 
we observe that, due to the last part of the proof of Corollary (TJ the approximate solution 
z t := [u e , ip e ] satisfies 

\\\u t (t)\\ 2 + v ( \\Vu e \\ 2 <]-\\uq\\ 2 + i l (aip t - J* ip e )u e ■ Vip e + [ (h,u e )dr. 
* Jo * Jo Jn Jo 

Passing to the limit as e — > in this last inequality and arguing as at the end of the proof 
of Corollary OQ taking advantage of (14742]) . ( fQgjl . (T47461) and ( 14747p . we get ( 13745]) . The 
energy inequality (I3.45P is then obtained by adding together (13.461) and (I3.48p . □ 
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It is worth noting that if the mobility degenerates slightly stronger, then the weak 
solution of Theorem [2] satisfies a physically more relevant energy inequality with respect 
to the energy identity (13.44j) . This is proven in the following 

Corollary 2. Let all assumptions of Theorem^ be satisfied and d = 2,3. In addition, 
assume that m'(l) = m'{— 1) = 0. Then the weak solution z = [u,<p] given by Theorem^ 
also fulfills the following integral inequality 



E(z{1)) + I (H|V„|| Z + 



J 



)dT<£(zo)+ / (h,u)dr. 



for all t > 0, where the mass flux J given by 

J = - m {cp)V{a(f - J*<p)- m(cp)F"(cp)Vcp 



(4.58) 



is such that 



J e L\Q T ) 



J 



€ L\Q T ) 



Proof. Setting J e := — y m e (ip e ) V/i e , from (14. 32ft we have ||J7^||i,2(q t ) < C and therefore 
there exists J G L 2 {Q T ) such that 



J € ^J weakly in L 2 (Q T ). 



(4.59) 



On the other hand, setting J e := — m e (</? e )V// e = ^/m e (ip e )J e , we have also || j7" e ||L 2 (Q T ) < 
C. Thus there exists J G L 2 {Q T ) such that 



J t -± J weakly in L 2 {Q T ). 



(4.60) 



Since m e (<p e ) — > m((p) strongly in U(Qt) for every r G [1, oo), from (14.591) we have 
\fm e ((pt)J e — ^ \/m((p)J weakly in L 2 "~ 7 (Qt) for every 7 G (0, 1] and therefore by com- 
parison with ( I4.60p we get 



J := y/m{ip)J. 

Let us observe now that inequality ( 14 . 2 9 [) can be rewritten in the form 

Je 2 



u\\ Vm, 



dT<£(z )+ I (h,u e )dr, Vi > 0, (4.61) 




y/m e ((p t 

where J e is defined in terms of <p e , that is, 

J e = -m e (v5 £ )V/i e = -m e ((p e )V (atp e - J * <p e ) - m e (yj £ )F"(yj e )Vy? 
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Using f H~46j) . (Oi| and f H~5?j) it is easy to see that 

J e ->> -m(^)V(a^> - J * <p) - m(^)F"(^)V<^ weakly in L 2 - r '(Q T ), Vr? G (0, 1). 

On account of (14.601) . we therefore see that J is given by 

J = -m(ip)V(atp - J * ip) - m(tp)F"(ip)Vip, (4.62) 

and that the last convergence is also in L 2 (Qt)- Hence, taking the inferior limit of the 
third term on the left hand side of (I4.6ip . we have 



J 



dr 



\jVdr < liminf 

e-S-0 



\J e \\ dr = liminf 



3c 



\fm t {ip e 



dr. 
(4.63) 



On the other hand, due to condition m'(±l) = 0, it can be shown (see [T3| Corollary, 
p. 417], cf. also [10, Theorem 2.3]) that the set {x G Q : |<£>(x,t)| = 1} has zero measure 
for almost any t G (0,T). Hence we have F e (<p e ) — > F(<p) almost everywhere in Qt, and 
by Fatou's lemma 



/ FUp)< liminf / F e (cp e 
Jn Ju 



(4.64) 



Finally, ( I4.58j) follows immediately from f)4.6ip on account of fl4.63p . f !4.64p . f !4.62p and 
the second weak convergence in (14 .42 p . □ 



Remark 6. In comparison with the analogous result for the case of constant mobility (see 
[T6| Theorem 1]), Theorem |2] does not require the condition |y? | < 1 (the assumptions 
on ip Q imply only that |^ | < 1). This is essentially due to the fact that here we are 
dealing with a different weak formulation. Therefore, if F is bounded (e.g. F is given by 
(I4.2p ) and at t — the fluid is in a pure phase, e.g., tpo — 1 almost everywhere in Q, and 
furthermore «o = u(0) is given in Gdiv, then we can immediately check that the couple 
[u, ip] given by 

u = u(-,t), ip = tp(-,t) = 1, a.e. in Q, a.a. t > 0, 

where u is solution of the Navier-Stokes equations with non-slip boundary condition, 
initial velocity field Uq and external force h, explicitly satisfies the weak formulation (14. 9p . 
(I4.10p . Hence, the nonlocal Cahn-Hilliard-Navier-Stokes model with degenerate mobility 
and bounded double-well potentials allows pure phases at any time t > 0. This possibility 
is excluded in the model with constant mobility since in such model the chemical potential 
\i (and hence F'(cp)) appears explicitly. 
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Remark 7. If m'(l) ^ and m'(— 1) ^ 0, then, as a consequence of (Al) and of the 
definition of the function M, both F and M are bounded in [—1,1]. This can be see, for 
instance, by writing M as follows 



In this case, the conditions F(ipo) G and M(<po) G L x (f2) of Theorem [2] are satisfied 

by every initial datum ip such that \(po\ < 1 in Q. Therefore the existence of pure phases 
is allowed. A relevant example for this case is given by ( 14. 2p . On the other hand, if 
m'(l) = m'{— 1) = (in this case we say that m is strongly degenerate), then it can be 
proved (cf. [T3], Corollary] and also [TQ1 Theorem 2.3]) that the conditions F(ipo) G L 1 (fi) 
and Af(^o) G imply that the sets {x G fi : v 9 o( a; ) — 1} an d {a; G f2 : ^(x) = — 1} 

have both measure zero. Hence, in this case we have obviously |<^ Q | < 1 and furthermore 
it can be seen that also the sets {x G Q : <p(x,t) = 1} and {x G Q : ip(x,t) = — 1} 
have both measure zero for a. a. t > 0. Therefore, in the case of strongly degenerating 
mobility the presence of pure phases is not allowed (even on subsets of Q of positive 
measure), and so, from this point of view the situation is more similar to the case of 
constant mobility. Summing up, as far as the possibility of existence of pure phases is 
concerned, the difference between the cases of constant and degenerate mobility is more 
relevant when the mobility is degenerate, but not strongly degenerate (on this issue see 



If, in addition to the assumptions of Theorem [21 we suppose that the initial chemical 
potential /x(0) =: /io belongs to H, then we can improve the regularity of \i. In particular, 
we obtain /i G L 2 (0,T; V). This is shown in the next theorem, which also requires some 
further reasonable conditions on the singular part F\ of the potential F (the potential (14.21) 
is however included) and whose proof is based on a suitable choice of the test function (see 
[22])- We point out that, while in [22] the improved regularity for the chemical potential is 
justified only by formal computations, the argument in the proof of Theorem [3] is rigorous. 

Theorem 3. Assume that all assumptions of Theorem^ are satisfied with J such that 



where k = 6/5 if d = 3 and k > 1 if d = 2. In addition, assume that Fx G C 3 (— 1, 1) and 
that there exist some constants a^, (3 > and p G [0, 1) such that the following conditions 
are fulfilled 





also gl]). 




pF['(s) + F%(s) + a(x) > 



Ws G (-1,1 



a.e. in f2 



(4.65) 



m(s)F['(s) > a , 



VsG [-1,1], 



(4.66) 
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\m 2 (s)F{"(s)\<p , Vsg[-1,1], (4.67) 

Fl(s)Fl"(s) > 0, Vse(-l,l). (4.68) 

Let fo be such that 

F'(<p ) G H. (4.69) 
Then, the weak solution z = [u,(p] given by Theorem^ fulfills 

fie L°°(0,T;H), V/i G L 2 (0, T; H). (4.70) 



As a consequence, z = [u,<p] also satisfies the weak formulation (13.51) and (I3.6p . the energy 
inequality (13. 7ft and, /or d = 2, the energy identity (I3.12p . 

Remark 8. Notice that in Theorem [3] assumptions (A2) and (A4) can be omitted. Indeed, 
(A2) follows from (14661) and (Al), while (A4) follows from fl4T65|) and fITOD . Notice 
moreover that the standard potential F and the mobility m in (14. 2p comply with the 
assumptions of Theorem [3J 

Proof. Let us write the weak formulation of (14.121) as 

(<p' e , if;) + (m e (</? e )V/T e , VV) + {m e (ip e )V(atp e + KM), W>) = (<We, VV), (4.71) 

for all if> G V, where we have set 

/T e := u> e + i^ e (^ e ), io e := - J * <£> e , 

and where F u is defined as in (I4.18I) . Take if) = F{ e ((p e )F" e ((p e ) as test function in (14.711) 
and notice that if) G V. Indeed, F{ e (ip e ) G V and, since by (14.181) F" £ is globally Lipschitz 
on M, we have also F" e ((p e ) G V. On the other hand, by means of an Alikakos' iteration 
argument as in [8, Theorem 2.1], we can see that (p e G L°°(Q). Hence F{ e ((p E ) and 
F£((p e ) are in V n L°°(fi). This yields if) = KM F uM eVn L°°(fi). Using this test 
function we first see that the contribution of the convective term vanishes. Indeed, by the 
incompressibility condition (11.81) we deduce 

/ (u e ■ V<Pe)Kfa)KM = [ «e ■ V(^|^) = 0. 

Furthermore, we can see that we have 

(^,KMKM) = \j t l \KM\ 2 - (4.72) 

Indeed, setting G e (s) := F[ 2 (s)/2, we find G" = F{ e F{" + F" 2 > almost everywhere in 
JR., due to (I4.68p . i.e., G e is convex on R. Hence, (I4.72p follows as an easy application of 
the chain rule result proved in |T2i Proposition 4.2]. 
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We are therefore led to the following identity 
1 d 



||F 1 > e )|| 2 + / m e (<f £ )F[' £ ^ £ )Vp e ■ VF' U {^) 
Jn 

+ / mMF{ € (<p e )F£(<p e )Vf}; e -V<p e + [ m e ( Ve ) Ve Va-V(F' u { Ve )F? e { Ve )) 
Jn Jn 

+ / m e (cp £ )(a + F&<p e ))V<p e ■ V(F{ £ (<p £ )F{' £ (<p £ )) = 0, (4.73) 
Jn 



which can be rewritten as 



\j t \\& ~ w t \\ 2 + J^m e {ip € )Fi' e {ip e )\V{p: e - w e )\ 2 + J m e (<p e )F? e {<p e )Vw e ■ V(/i e - w e ) 

+ / m £ (ip e )(jI € -w e )Fl"(ip e )VjIe-Vipe+ / m e (ip e )ip e F" e (ip e )Va-VF[ e (ip e ) 
Jn Jn 

+ / m e (<p e )<p e F' u (<p e )F£(<p e )Va-V<p e + [ m £ {ip £ )(a + KM)Fi? (^)|V^ e | 2 
Jn Jn 

+ [ m £ (ip £ )(a + F^ £ ))F[ £ {<p £ )F^ £ )\V^ = 0. (4.74) 
Jn 

Condition (I4.65P implies that 

pF[' £ (s) + F£(s) + a(x) > 0, VsGl, a.a i6fi, Ve G (0, 1). (4.75) 
Then we get 



m e {y e )F? £ {y £ )\V{&-w e )\ 2 + / me(<p e ){a + FL(<p £ ))F£(<p e )\V<p e 
n Jn 



2 



> (1 - p) I m e (^)^' e (^)|V(/Ie - w)\ 2 . (4.76) 
n 



On the other hand, condition (I4.68P implies that, for all e G (0, 1), F[ £ F" £ > on R. Thus 
using again (14 .75 p . we have 

/ m e ^ e )(jI £ -w £ )F{'^ e )Vfi: e .V Ve + I m e ( Ve )(a + F^ e ))F{ £ ( Ve )F{'^ e )\V^ e \ 2 
Jn Jn 

> [ m £ {if £ )F[ £ ^ £ )F^ e )Vp £ -V^- p [ m e (^)Fr e (^ e )^ e (^X(^)|V^| 2 
Jn Jn 

m e {^)F[ e {^)F'l' £ { Ve )F^ £ )\V^ + [ m e (<f e )Fi £ (<f e )F^(ip e )Vw e • 
n Jn 



- P / m e (^)F^jF^F^|V^| 2 

= (l-p) [ m £ (<p £ )F[ £ ^ £ )F^ e )V(p £ -w £ )-V^ e + [ m&jFlMFgfa) Vw £ ■ Vy e 
Jn Jn 

= (l-p) / m e (^ e )F[ e (^ e )F[^ e )V/r e .V^ + p / m e (^ e )F[ e (^ e )^(^ e )V«; e • V^ e . 



(4.77) 
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Furthermore, observe that 



where 



* W : = ™.(.«(.)' % € 
Therefore the following identity holds 

/ m e (<^V e F^ e )F£V)Va • Vip e = [ j e (ip e )(p; e -w e )<p e Va-V(fi; e -w e ). (4.79) 
Jn Jn 

By plugging (14. 76 p . (I4.77P and (14.791) into (I4.74p and using ( I4.78P once more, we get 



(4.80) 



2^11^ - Well 2 + i 1 ~ P) J^ m ^e) F "e(Ve)\y(He ~ W £ )\ 2 

+ / m e (ip e )F" e ((p e )Vw e • V(/T e - W e ) + (1 - P) / J e (Ve)(fie ~ W e )\V(jX e - W £ ] 

Jn Jn 

+ / 7 e (v3e)(/ie - w e )V{fI e - w e ) • Vw e + / m e (^ e )F" e (^ e )</? e Va • V(/T e - w e ) 
Jn Jn 

+ / %(<p e )(p: e - w e )(p e va -v(p; e - w e ) < o. 

Jn 

Notice that the fourth term on the left hand side is nonnegative since there holds 

On the other hand, condition (14.661) implies that 

m e (s)F{' e (s) > oq, Vs 6 1, Ve G (0, 1), (4.81) 
and, using (H2), we have 

/ 7e(^)(/Ie - W e )V(jI e ~ W e ) ■ Vw e > -(1 - p)^||V(/T e - W e )f 

-d I 7e 2 (^)|V We | 2 |/I e - W£ | 2 , (4.82) 



7 e (V3 e )(/i e - W e )(p t VCL ■ V(/4 - W e ) > -(1 - p) — ||V(/i e - W e 

o 



-d / 7e 2 (^)^|Va| 2 |/Ie-We| 2 , (4.83) 
Jn 

/ m e (tp e )Fi' e (<pe)<feVa-V(p; e -w e ) > -(1 - p)^|| V(/T e - w e )|| 2 - C 2 ||^|| 2 , (4.84) 

Jn o 
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/ m e (<p e )F&(cp e )Vw t ■ V(/T e - w e ) > -(1 - p)^|| V(JT e - w e ) f - C 3 ||V™ e || 2 . (4.85) 
Jn 8 

In (I4.84p and (I4.85P we have used the fact that (14.531) holds as a consequence of the 
boundedness of mF" on [—1, 1]. 

Recall that %(s) = for all \s\ > 1 - e, and 7 e (s) = m 2 (s)F"'(s)/m(s)F"(s) for 
\s\ < 1 — e. Hence, on account of (I4.66P and (14.671) . we have 

|7e(s)|<7o, VsGR, Vee(0,l), (4.86) 

with 70 := Pq/(Xq. Thus the second term on the right hand side of (14.831) can be estimated 
as follows 

/ 7 e W'l Va| 2 |/I - w e \ 2 < 7o 2 || Vo||ioo (n) ||/z e - w £ \\ 2 . 

J{\<p e \<l-e} 

By means of (I4.8ip - (l4.87p . from (14.801) we deduce the inequality 

d 



(4.87) 



— N/ x e - w e || 2 + a (l - p)||V(/i e - w £ )|| 2 

< c 4 ||v^ e || 2 + c 5 / {i + \vw e \ 2 ) \jr e 

Jn 



w e \ 2 -w e \\ 2 + Cell^f. (4.t 



Observe now that 



C 4 / \Vw e \ 2 \fi e -w e \ 2 <C 4 \\iJ e -w e \\ L 6 {n )\\Vw e 
Jn 

< y (1 " P)ll V(/T e - w t )\\ 2 + C 7 (l + ||V We || 4 i6(n) ) ||£ - w t \\ 2 , (4.89) 

where we have supposed that d = 3. If d = 2, (14.891) still holds with the norm || Vtu e ||.£,6(Q) 
replaced by ||Viw e ||ir(n) for any given r > 4. 

On the other hand, by Young's inequality we have 

||Vu> e || L 6 (Q) < iV 3 ||^ e ||£,8(n), d = 3, 

||VW e || L 4 K /(4-3K) (n) < iV 2 ||^ e ||jr,4(0), d = 2. 

Since {y>J is bounded in L 2 (0, T; L 6 (fi)) n L°°(0, T; L 2 (ft)) L 4 (0, T; L 3 (fi)) if d = 3, 
and bounded in L 2 (0, T; V) n L°°(0, T; L 2 (ft)) <-> L 4 (0, T; L 4 (fi)) if d = 2, independently 
of e (cf. (I4.6P ). then, using also (H2), we get 

J 1 1 11x6(^)^1" < C, d = 3; J \\Vw e \\U/(4- 3KHn) dT <C, d = 2. (4.90) 

Combining (I4.88P with (I4.89p . and using also (14.311) . we are led to the inequality 

j t \\& ~ ™ e || 2 + y (1 " P)I|V(£ - w e )\\ 2 < a + ,j e ||£ - W J 2 , (4.91) 
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where a is a positive constant and where the function r) e , given by r\ t := 6*7(1 + 1| Viy e ||^ 6 (n)) 
if d = 3, or by r) e := Cs(l + || Viy e ||^4 K/( 4_ 3K) ^) if d = 2, on account of (I4.90p . satisfies 

It is easy to see that 

K00I < 1^001, Vs G (-1,1), Ve G (0,e ], (4.92) 
with €q > small enough. Therefore, condition (I4.69P entails 

||/r 6 (o)-w 6 (o)|| = 11^0)11 < WKMW, 

and by applying Gronwall's lemma to (I4.9ip and using the fact that w e is bounded inde- 
pendently of e in L°°(0, T; H) as well as its gradient (cf. (H2) and (|4.6j) ). we immediately 

get 

fie L°°(0,T;H), VJle L 2 (0,T;H), (4.93) 
where Jl := w + Fi(ip) and w := — J * (p. Furthermore, due to (14.191) . we have (see also 

\F^(s)\ < L 3 \s\+L 4 , \Fg(s)\ < L 5 , Ws G R, 

for some nonnegative constants L 3 , L 4 , L 5 which are independent of e. By using these 
estimates, (14.3 ip and (I4.35P we also obtain 

F* 2 {<p) G L°°(0, T; H), VF^ip) G L 2 (0, T; H). (4.94) 

From (1433]) and (fl~53|) we immediately get (14770]) . 

In order to see that (I3.5P and (I3.6P are satisfied, consider the weak formulation (I4.5ip . 
(I4.52p of the approximate problem P e , write F"((p e )'Vip e = Vi^'^e), rearrange the terms 
on the left hand side of (14.511) . take the definition of \i t into account and pass to the limit 
as e — > 0. 

The energy inequality (13. 7p can be deduced by passing to the limit in (I4.29P and by 
observing that 

/ F e (<p ) < [ F^cpo) + [ F 2e {cp ), 
Jn Jn Jn 

and that F 2)E (<y2o) — > J n F 2 (v9 ) by Lebesgue's theorem (see (I4.23P ). 

Finally, the energy identity ( 13 . 1 2[) for d = 2 is obtained by choosing if) = \i as test 

function in (13.51) . □ 
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The argument of Theorem [3] can be exploited to prove rigorously the existence of a 
weak solution satisfying (14. lip with 6 < p < oo, for d = 3. We point out that the L°°(L P ) 
regularity of (p follows only formally from the condition \(p\ < 1 proved in Theorem [21 
Indeed, in the case d — 3 it is not known that the if component of the weak solution 
of Theorem |2j is measurable with values in LP if p > 6. Instead, a rigorous way to 
deduce such regularity is by means of an approximation argument which makes use of an 
approximating potential with polynomial growth of order p. 

Corollary 3. Let all the assumptions of Theorem^ and Theorem^ be satisfied and d = 3. 
Ln particular let ipo be such that (I4.69P holds. Ln addition, assume that there exists €q > 
such that the following assumptions are satisfied 

F"'(s) > 1, Vs G [1 - e , 1) and F?'(s) < -1, Vs G (-1, -1 + e ], (4.95) 
F[(s) > 0, Vs G [1 — e , 1) and F{(s) < 0, Vs G (-1, -1 + e ]. (4.96) 

Let p G (6, oo) be fixed arbitrarily. Then, for every T > 0, there exists a weak solution z = 
[u,<p] to ( I1.5p - (jl.l0p on [0, T] in the sense of Definition^ satisfying the weak formulation 
( I3.5p . (13. 6p and such that 

<p G L°°(0, T; L P (Q)). (4.97) 



Proof. We argue as in the proofs of Theorem |2j and Theorem [31 but we employ a different 
approximation of the singular part F le of the double-well potential. Namely, here we 
consider the approximate problem (I4.12p - (l4.17p with the following choice for F le 

,2 



Fu{s) 



[ F 1 {l-e) + F[(l -e)(s-(l- e)) + fif (1 - e)(s - (1 - e)f 
+l(s-(l-e))\ s>l-e 
Fi(s), |s|<l-e 

F x (-1 + e) + F{(-1 + e)(s- (-1 + e)) + \F^-l + e)(s- (-1 + e )) : 



(-l + e) 



s < -1 + e. 



Let us argue as in Theorem [2] to deduce the estimates for ip e and u e (and their time 
derivatives). Instead, the passage to the limit will be performed differently with respect 
to Theorem [2j 

We first check that ( I4.2ip is still true for this new approximate potential. Indeed, for 
1 — e < s < 1, thanks to (I4.95P we have, for some £ G (1 — e, s) 

F 1 (s) = F 1 (l -e) + F{(1 - e)(s - (1 - c)) + - e)(s - (1 - e)) 2 
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>F 1 (l-e) + F[(l -e)(s-(l- e)) + - e)(s - (1 - e)) 2 

+ ^(^-(l-e)) P = Fi,(3), 

provided that e G (0, e ], with e as in (14.951) . For — l<s<— 1 + ewe argue similarly. 

It is easy to see that the uniform (with respect to e) coercivity condition (I4.25P will 
now be replaced by (cf. also (A2)) 

F e {s) > ^\s\ p -c p , VsGE, Vee(0,eo], (4.98) 

where the constant c p depends on p but is independent of e. It is also immediate to check 
that (I4.28P and (H7) hold for the new approximate potential. Therefore, since we have 
also F e G C^ C (R), assumptions (H1)-(H6) of Theorem [TJ with (H4) replaced by (H7), are 
satisfied. Hence, for every fixed e G (0,eo], Problem P e admits a weak solution satisfying 
(133]) . the first of (JHZ1, (ED, (13TT0]) . and the energy inequality (jOg) . By means of 
(I4.28P and (14 . 9 8 [) . from the approximate energy inequality we can still recover estimates 
fl430|) . fl432|) . ( jQ5D and fl436|) . while estimate (Olj) will now be substituted with the 
stronger 

||^e|U°°(0,T;iP(n)) < C. (4.99) 

The estimates for u' e and ^ can be improved as well. Indeed, instead of (I4.38P and (14.411) 
we now find 

W e \\&(aW)<C> d = ^ ( 4 - 10 °) 

and 

KIU 2 (o,T;VlJ < C> II^IU^o^v") < C, d = 2, (4.101) 

respectively, while (I4.4UI) still holds. By compactness results we can still deduce that there 
exists a couple z := [u,ip] and a subsequence such that (I4.42I) - (I4.49I) hold. Instead, in 
(g35D, (g3HD and fl4l9l we shall have L 2 in place of L 2 " 7 , L 4 / 3 and L 2_<5 , respectively, 
and 

<p e ->> <p weakly* in L°°(0, T; L p (fi)), weakly in L 2 (0, T; V), (4.102) 
in place of fT4~46]) . 

As far as the passage to the limit is concerned, the situation is now a bit differ- 
ent. Indeed, due to the polynomial growth of Fu, a strong convergence for the term 
m e (ip e )F{' e ((p t ) in some L s (Qt) space with s > 2 is no longer available in order to de- 
duce the weak formulation (14. 9p for (p. Therefore, the idea is to use (I4.69p . argue as in 
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Theorem [3], get a control for /x £ in L 2 (0,T; V) and pass to the limit to deduce the weak 
formulation ( 13. 5 \ for (p. Of course, we have to check that the argument of Theorem [3] still 
applies with this different choice of Now, it is easy to verify that the test function 
if) = Fl e ((p e )F" e ((p e ) in the weak formulation of (I4.12p is still in V. Indeed, F" e is now 
locally Lipschitz on R and tp e E V n L°°(Q), which implies that F" € ((p e ) E V n L°°(Q). 
Thus ip EV. 

By repeating the calculations in the first part of the proof of Theorem [3] we are then 
led again to the differential identity ( I4.74p . It is immediate to check that ( I4.8ip still holds. 
As far as (I4.86P is concerned, notice that in this case we have for, e.g., s > 1 — e 

% S rn{\ - e) (if (1 - e) + Jp(p - 1) (a - (1 - e))^ 2 ) 

< Cl|P m 2 (l-e)( a -(l-e)) p - 3 
«o + c 2 , p m(l - e) (s - (1 - e)) P 2 

where (I4.66P has been taken into account. Now, if we consider the function 

S 2 p-3 

i 

by an elementary calculation we can see that ips has its maximum at r* := ((p — 3)/5) p_2 
and that ^(r) < ^(r*) = c^fr- 1 ^ -2 ) < c p , for all r > and for all 5 E (0, 1), where 
c p is a positive constant depending only on p. Hence we deduce that |7 e (s)| < c' for all 
\s\ > 1 — e and for all e E (0, 1). Using also the fact that %(s) = m 2 (s)F"'(s)/m(s)F"(s) 
for \s\ < 1 — e, and taking (I4.66P and (14.671) into account, we can therefore conclude that 

\le{s)\ < 7o, P , Vs £1, Ve E (0, 1), (4.103) 

with 7o iP given by 7 0iP = max(c p , 0o/ a o) an d independent of e. Hence, (I4.86P is still 
satisfied. 

On the other hand, using (A2) and (I4.96p . it is easy to check that also conditions 

KeK > on R, Ve E (0,1), 

and (I4.92p still hold. We can therefore conclude that the argument of Theorem [3] applies 
and that, due also to (I4.69p . the following estimates hold (cf. (I4.9ip ) 

||^e|U«(0,r;fl) < C, ||V/i e ||L2(0,T;H) < C, (4.104) 

where the constant C depends on \\fio\\ or, equivalently, on ||F'((^o)||- 

Finally, let us consider the weak formulation (I3.5P and (I3.6P for the approximate 
problem (I4.12p - (l4.17p . By passing to the limit as e — > and by means of (14.1041) and the 
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convergences provided in the proof of Theorem [2], it is not difficult to show that z = [u, (p], 
which in particular satisfies ( I4.97P due to the first of (14.1021) . satisfies the weak formulation 
dS3D and (USD. □ 



5 The global attractor in 2D 

In this section we consider system fll.5j) — fll.lQp for d = 2 and we suppose that the external 
force h is time-independent, namely, 

(A5) h g V' m . 

Let us introduce the set Q mo of all weak solutions to (I1.5I) - (I1.10I) (in the sense of 
Definition [2]) corresponding to all initial data z = [uo,(fo] G X mo , where the phase space 
X mo is the metric space defined by 

"^roo • ^ div X y mo > 

with y mo given by 

y mo ■= {<p e L°°(tt) : \<p\ < 1 a.e. in ft, F(ip) e L\n), \lp\ < m }, (5.1) 

and m G [0, 1] is fixed. The metric on X mo is 

d(z 2 ,zx) := \\u 2 - Ui\\ + ||v?2 - <fi\\, 

for every z x := [ui,v?i] and z 2 ■= K,^] in X mo . 

In the next proposition we shall prove that the set of all weak solution is a generalized 
semiflow in the sense of J.M. Ball (cf. [7]). 

Proposition 1. Let d = 2 and suppose that the conditions of Theorem^ are satisfied. 
Assume also that (14.651) . (I4.66P and that (A5), (H2) hold. Then Q mo is a generalized 
semiflow on X mo . 

Proof. It is immediate to check that hypotheses (HI), (H2) and (H3) of the definition 
of generalized semiflow [3 Definition 2.1] are satisfied. It remains to prove the upper 
semicontinuity with respect to initial data, i.e., that Q mQ satisfies (H4) of Definition 
2.1]. Take then Zj = [uj, ipj) G Q mo such that Zj(0) — > z in X mo . Our aim is to prove that 
there exist z G Q mo with z(0) = z and a subsequence {zj k } such that Zj k (t) —> z(t) in 
X mo for all t > 0. Now, each weak solution Zj = [uj, <pj] satisfies the regularity properties 
fO]) -(ra. $4~TT\ and the energy equation fl3T44"j) which, on account of fl465|) . (jMSD and 
of the fact that \tpj\ < 1, implies 

j t (INI 2 + Itoll 2 ) + (i - p)«oI|v^-|| 2 + H|v % || 2 < c + c \\ Uj r + ^\\Hv L , 
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where the positive constant c depends on J and on m. By integrating this inequality 
between and t, using the fact that Zj(0) — > z$ in X mo and Gronwall's lemma, from the 
above differential inequality we get 

\\Uj\\L°°iO,T-,G div )nL2(0,T;V div ) < C, \\ (fj \\l™(0,T;H)iiL2(0,T;V) < C. 

By comparison in the variational formulation (14. 9 p and ( I4.10p written for each weak 
solution Zj = [uj, (fj] we also obtain the estimates for the time derivatives u'a and <p'j 

\\u'j\\ L Ho,T;V^ v ) < C, \\w'j\\L2{o,T-y') < C. 

Therefore, by standard compactness results, we deduce that there exist two functions u 
and ip such that, for a not relabeled subsequence, we have 

Uj — u weakly* in L°°(0, T; G div ), weakly in L 2 (0, T; V div ), (5.2) 
Uj w strongly in L 2 (0, T; Gdiv), a.e. in Q T , (5.3) 



■it 



u t weakly in L 2 (0,T;^J, (5.4) 

^ ->> <p weakly* in L°°(0, T; weakly in L 2 (0, T; V), (5.5) 

<fj^-<fi strongly in L 2 (0, T; if), a.e. in Q T , (5.6) 

V?j ->■ (ft weakly in L 2 (0, T; V). (5.7) 

Passing to the limit in the variational formulation for Zj = [uj, tpj] and using the weak/ strong 
convergences (I5.2p - (l5.7p and the fact that \<pj\ < 1, we get that z := [u, ip] G Q mo (the 
argument is similar to the passage to the limit in (I4.5ip and (I4.52p in the proof of Theorem 
|2]). Furthermore we have z(0) = z Q , since, as a consequence of ( 15. 2p . ( 15. 4p . ( 15. 5 p and ( 15 .7p . 
we have, for all t > 

Uj(t) — ^ u(t) weakly in Gdiv, (5.8) 
(Pj(t) — ^ (p(t) weakly in H. (5-9) 

Let us now see that Zj(t) — > z(t) in X mo for all t > 0. First observe that the energy 
equation (13.441) can be written in the form 



d_ 

di 



where 



E (zj(t)) + / rn((pj)F"(<p j )\V<p j \ 2 + / am(^)|V^f + ^||V^|| 2 = 0, (5.10) 
Jn Jn 

E( Zj (t)) :=E( Zj (t)) - [ [ m(^)(VJ*^-^Va)-V<^ 

Jo Jn 

— / (aifj — J * <Pj)uj ■ Vfj — I (h,Uj)dr, 
Jo Jn Jo 
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and 

E{ Zj {t)) :=\(\\um\ 2 +\wm\ 2 )- 

On the other hand, from (15. 3p and (15.61) we have, for almost any t > 0, 

Uj(t) -)■ strongly in 

(fij(t) —> (f(t) strongly in H, 

and hence E(zj(t)) — > E(z(t)) for almost any t > 0. Moreover, due to the condition 
1^1 < 1 for each j, to the pointwise convergence (15.61) and Lebesgue's theorem we have 

ipj — > if, strongly in L s (Qt), Vs G [2, oo). 

Therefore, recalling that m(ipj) — > m(ip) strongly in L s (Qt) for all s G [2, oo) and using 
the second weak convergence in (15.51) . we have, for all t > 0, 

/ / m(ipj)(yj * ifij — (fjVa) -Vifij — > / / m(ip) (VJ * <p — ipVa) ■ Vy2. 
Jo Jn Jo Jn 

Furthermore, since a<pj — J * ipj a<p — J * <p strongly in L s (Qt) for all s G [2, oo) and, 
by (15. 2p . uj — > u strongly in L S {Q T ) 2 for all s G [2, 4), then we also have, for all t > 

/ / (aifj — J * <fj)uj ■ — > / / (aip — J * (p)u ■ V(p. 

Jo Jn Jo Jn 

Therefore 

E(zj(t)) E(z(t)), a.a. t > 0. 

Since, due to (15.101) . E(zj(-)) is non increasing in [0, oo) for every j and £ , (^(-)), E(z(-)) G 
C([0, oo)), we infer 

E{zj{t)) -+E(z(t)), Vt>0, 

and hence also E(zj(t)) — > E(z(t)) for all t > 0. This last convergence together with 
(I5.8P and (15.91) yield Uj(t) — > u(t) strongly in G^w and <fj(t) — >■ <^(t) strongly in if, and 
hence z,(t) — > z(t) in X mo , for all t > 0. □ 

Proposition 2. Lei d = 2 and suppose that the conditions of Theorem^ are satisfied. As- 
sume also that (14.651) . (I4.66P and (A5) hold. Then Q mo is point dissipative and eventually 
bounded. 
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Proof. Let us estimate the second term on the right hand side of (13.441) as follows 
/ {cap — J * <f)u ■ Vp> = I (-ip 2 Va — (p(y J * ip)) ■ u 

<^|Q|V2 M <^|| Vw ||2 + C(o) 

where Co = 95 2 |f2|/4i/Ai with the constant b defined as in (H2). Moreover, the first term 
on the right hand side of (13.441) can be controlled in the following way 

/ m(ip)(VJ*p-<pVa) • V<p < 2m |ft| 1/2 6|| V<p\\ < (1 - p)^|| V^|| 2 + d, 
Jq 2 

where C\ = 2ml\Q\b 2 / cxo, and m* = max sg [_ ljl ] m(s). Then, by taking (I4.65[) and f )4.66[) 
into account, we get the differential inequality 



d 



u 



+ |M| 2 ) + (1 - p)a ||V^|| 2 + is\\Vu\\ 2 < C 2 + 



km, , 

7/ div 



where C*2 = 2(Cq + C\). By using the identity ||<p|| = \\ip — <p \\ + \\(p \\ , we obtain 



d 
~d~t 



u 



Iv 5 - ^of) + (1 - p)®oCp\\<p - <Pq\\ 2 + ^Ai||m|| 2 < C 2 + -\\ h \\v> j 



where Cp is the constant appearing in the Poincare-Wirtinger inequality. Therefore we 
get 

d 



dt 



E(z(t))+ V E(z(t))<C 3 , 



(5.11) 



where z = [u, ip— <p> ] and the constants r\ and C3 are given by r\ = min((l— p)aoCp, u\i)/2, 
and C3 = C2/2 + \\h\\y, /2u, respectively. By Gronwall's lemma from (15.1 ip we have 

div 

E(z{t)) < E(z(0))e-* + — , Vt > 0. 
This estimate easily yields 

2C 

d 2 (z(t), 0) < d 2 (z , 0)e-"* + — + |^ | 2 |f2|, Vt > 0, 

V 

which entails both the point dissipativity and the eventual boundedness of Q mo - 
We can now state the result on the existence of the global attractor. 



□ 



Proposition 3. Let the assumptions of Proposition^ hold. Then Q mo possesses a global 
attractor. 
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Proof. Thanks to Proposition |2] and by Proposition 3.2] and [7J Theorem 3.3] we only 
need to show that Q mo is compact. Let {zj} C Q be a sequence with {%(0)} bounded in 
A" mo . We claim that there exists a subsequence {zj k } such that Zj k (t) converges in X mo 
for every t > 0. Indeed, the energy equation (13.44)) entails the existence of a subsequence 
(not relabeled) such that (cf. the proof of Proposition [1]) , for almost all t > 0, 



(t) — > u(t) strongly in Gdiv, ~^ vify strongly in H and a.e. in Q, 



where z = [u,tp] is a weak solution. By arguing as in the proof of Proposition [TJ we infer 
that E(zj(t)) E(z(t)) for all t > 0. Thus Zj(t) z(t) in X mo for t > 0, which yields 
the compactness of Q mo . □ 

Remark 9. We point out that the existence of the global attractor has been established 
without the restriction |^| < 1 on the generalized semiflow (compare with Remark [6]). In 
particular, this result does not require the separation property. It is also worth observ- 
ing that in dimension three the trajectory attractor approach used in [T5J [16] might be 
extended to the present case. 

6 The convective nonlocal Cahn-Hilliard equation 
with degenerate mobility 

By relying on the results of the previous sections we can prove similar results for the 
convective nonlocal Cahn-Hilliard equation with degenerate mobility and with a given 
velocity field, for d — 2, 3. In particular, from Theorem [2] it is straightforward to deduce 
the following 

Theorem 4. Assume that (Al)-(M) and ( H2 ) hold. Let u G Lf oc ([0, oo); Vdiv H L°°({l) d ) 
be given and let (po e L°°(f2) be such that F(ip Q ) G L 1 (fi) and M(ip ) G L l (Q). Then, for 
every T > 0, there exists a weak solution <p to fll.5p . (11. 6p . (I1.9p 1 on [0, T] corresponding 
to ifo fulfilling ( I4.6p ( 1478]) . the weak formulation ( 14 . 9 )) . and such that ^{t) = Tp^ for all 
t G [0, T\. Furthermore, p satisfies 

tpeL°°(0,T;L*>(n)), (6.1) 

where p < 6 for d = 3 and 2 < p < oo for d = 2. In addition, for almost any t > 0, the 
following energy identity holds 

Id 2 
2dt M + > m 



((^)F"(^)|V^| 2 + / am(ip)\V<p\ 2 + / m(y?) (y?Va - V J * ip) ■ Vp = 0. 
Jn Jn 

(6.2) 
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Remark 10. The basic estimates in the proof of Theorem H] are now obtained from (16. 2p 
(written for the approximate solution cp e and integrated between and t), in place of 
( gjZgp . In particular, (JO} still yields ( lOTj) . (Q5|) and (OB) , so that we can argue as 
in the proof of Theorem [3J 

For the convective nonlocal Cahn-Hilliard equation with degenerate mobility we can 
also prove a uniqueness result (see [21] for more restrictive assumptions). We remind that 
uniqueness of solutions is an open issue for the local case (see [13]) as well as for the 
complete system fll.5p - fll.10p even in dimension two. 

Proposition 4. Let all the assumptions of Theorem^ and ( I4.65p . (I4.66P be satisfied. 
Then, the weak solution to ( II. 5p . fll.6p . f 1 1 . 9 p i . (ll.lOp ? is unique. 



Proof. Following [21], let us introduce 



Ai(s) := / m(a)F['(a)da 7 A 2 (s) := / m(a)F^'(a)da, T(s) := / 
Jo Jo Jo 



m(a)da, 



for all s e [-1, 1]. Due to (Al) and (ITOj) we have Ax e ^([-1, 1]) and < a < A[(s) < 
a for some positive constant a. Then, it is easy to see that the weak formulation (14.91) 
can be rewritten as follows 

(ip t ,i/;) + (VA(-,<f), W) - (r(^)Va, W) + {m(<p)(<pVa - VJ * <p), Vtp) = (u</?,Vf), 

(6.3) 

for all ip G V, where 

AO, s) : = Ai(s) + A 2 (s) + a(x)T(s). 

Consider now (16. 3p for two weak solutions tp\ and </2 2 . Let us take the difference between 
the two identities, set ip := ipi — <y2 2 and choose ip = Afip as test function in the resulting 
identity (notice that Tp = 0). This yields 

\j t \W V M 2 + (Afo) - Afa), v) - (W - %))Va, V/V <p) 

+ ((m(</? 2 ) — "^(^OXv^Va — VJ * y? 2 ) + m(y>i)(</?Va — VJ * <p), V/Vy?) 

= (u^,VA/>). (6.4) 

On account of (I4.65P and (I4.66p . we find 

(A(-, <p 2 ) - A(-, (p) > (1 - p) / m(^ 2 + (1 - 6)^)^(6^ + (1 - #)y?i V 2 

./n 

> (l-p)a |M| 2 . (6.5) 
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Furthermore, since \<fii\, \(p%\ < 1, then we have 

| ((Tfoa) - r(^))Va, VA/» | < ||V^|| < ~(1 - p)a |M| 2 + C x \\VN^pf, 

(6.6) 

\{m(tp 2 ) -m{ip 1 )){ip 2 Va-VJ*ip 2 ),VN<p) \ < 2m„6||y?|| || VA/>|| 
<i(l-pH||^|| 2 + C 2 ||V^|| 2 , (6.7) 
| (m(v?i)(^Va - VJ * <p), VNip) | < 2m*fo||^|| \\VJ\f<p\\ 

< hi - p)a y\\ 2 + C 3 \\V^\\ 2 , (6.8) 



1 

lL°°(Q) d 



< -(1 - p)a |M| 2 + ^llttlHoofn^H^/Vir, (6-9) 



where b is given as in (H2), m, as in the proof of Proposition m** := max se [_! ^ m'(s) 
and the positive constants C\,..,C± depend on a 0) P, m *i m **- By plugging (I6.5p - fl6.9p 
into ( 16. 4p we get 

j t \W V M 2 + (1 - p)aoh\\ 2 < C 5 (l + N||. (n) ,) II^VII 2 (6-10) 
and Gronwall's lemma applies. □ 



As a consequence of Theorem H] and of Proposition H] we can define a semiflow S(t) on 
3^m ( c f- (15. ip )- mo G [0,1], endowed with the metric induced by the L 2 — norm. It is then 
immediate to check that the arguments used in the proofs of Proposition [TJ Proposition 
[2] and Proposition [5] can be adapted to the present situation. This gives 

Theorem 5. Let (A1)-(A4), fl4~65j) and fl4~66]) hold. Suppose that u G L°°(Q) d is given 
and independent of time. Then the dynamical system (y mo ,S(t)^ possesses a connected 
global attractor. 

Remark 11. The fact u must be divergence-free is not necessary (see also [751 [TB"] ). In 
particular, the convective term can be of the form V • (u<p). 

Remark 12. We stress once again that this is, to the best of our knowledge, the first 
result on the existence of the global attractor for the nonlocal Cahn-Hilliard equation 
with degenerate mobility and singular potential. 
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